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transition. 
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I. INTRODUCTION 

Over the past decade holography has emerged as a valuable tool for gaining insight into the 
physics of strongly coupled gauge theories. Top-down studies based on string/M-theory setups 
have been met by a number of bottom-up constructions, with applications ranging from the 
realm of quantum chromodynamics (QCD) to that of condensed matter systems. Holographic 
techniques have been particularly useful for probing the hydrodynamic regime of strongly in- 
teracting thermal field theories, a regime which is notoriously difficult to study directly and - 
unlike thermodynamics - poses a challenge to lattice simulations. 

Within this program, many efforts have been directed at better understanding the dynamics 
of the strongly coupled QCD quark-gluon plasma (QGP), and in particular at computing its 
transport coefficients. One of the most exciting results which has emerged from the heavy ion 
program at RHIC - and now at LHC - is the observation that the hot and dense nuclear matter 
produced in the experiments displays collective motion. In fact, the QGP fireball created in off- 
central collisions is not azimuthally symmetric, but rather shaped like an ellipse. As a result, the 
pressure gradients between its center and its edges vary with angle, giving rise to an anisotropic 
particle distribution. The matter formed in the collisions then responds as a strongly coupled 
fluid to the differences in these pressure gradients, displaying a collective flow which is well 
described by nearly ideal hydrodynamics with a very small ratio of shear viscosity to entropy 
density ^. 

Experimentally, the flow pattern can be quantified by Fourier decomposing the particles' 
angular distribution. In particular, it is the second Fourier component V2, the so-called elliptic 
flow, which is the largest in non-central collisions and is the observable most directly tied to 
the shear viscosity. Thus, bounds on ^ can be extracted from elliptic flow measurements, with 
the most advanced analysis at the moment giving 4ir 2 < 2.5 [1]. Higher order harmonics 
- initially neglected because they were assumed to be too small for symmetry reasons - also 
play an important role in determining the shear viscosity (see Section IV for a more detailed 
discussion). We refer the reader to [2, 3] for some early references on the RHIC results and the 
range of -, and to [4-7] for more recent ones including discussions of the first LHC results. 
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A remarkable result that has emerged from holographic studies of strongly coupled gauge 
theories has been the universality of the shear viscosity to entropy ratio [8, 9], which was shown 
to take on the particularly simple form ^ = 4- in any gauge theory plasma with an Einstein 
gravity dual description 1 . Its order of magnitude agreement with RHIC (and now LHC) data 
was one of the driving motivations behind the efforts to apply holography to the transport 
properties of the QGP (see [11, 12] for recent reviews). It is by now well understood that 
deviations from the universal result 2 = -L (both below and above) are generic once curvature 
corrections to the leading Einstein action are included 2 . Moreover, when there is another scale 
A in the system in addition to temperature T, the viscosity to entropy ratio typically runs as 
a function of T/A in such higher derivative theories. We should emphasize that this type of 
temperature flow for the shear viscosity arises in a number of holographic constructions, from 
theories of higher derivatives in the presence of a chemical potential [16-18] or non-trivial scalar 
field profiles 3 [20, 21] and also to systems with spatial anisotropy [10, 22-27]. 

The viscosity to entropy ratio is in fact known to be temperature dependent for a variety 
of liquids and gases in nature (as well as for ultracold fermionic systems close to the unitarity 
bound), exhibiting a minimum in the vicinity of a phase transition (see Figure 1). A similar 
behavior is expected [28] for 2 near the temperature T = T c of the QCD phase transition which 
separates hadronic matter from the QGP phase. In the hadronic phase below T c , the hadronic 
cross section decreases as the temperature is lowered, leading to an increase in ^ [29, 30] (for an 
analysis of transport in the hadronic phase see e.g. [31]). On the other hand, in the deconfined 
phase at temperatures well above T c , asymptotic freedom dictates that 2 should increase with 
temperature (the coupling between quarks and gluons decreases logarithmically [32, 33]). From 
the behavior in these two opposite regimes, we conclude that we should expect a minimum for 
- somewhere in the intermediate range. 

s ° 

A precise determination of the temperature behavior of transport coefficients such as 2 is an 
important ingredient for understanding the dynamics of the strongly coupled medium produced 

1 An exception is the case of anisotropic fluids, as first observed in [10]. 

2 For reviews of the shear viscosity bound we refer the reader to [13-15]. 

3 Dilatonic couplings to higher derivative terms in the context of ^ have also been studied in [19] . 
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FIG. 1: Schematic plots of the shear viscosity to entropy density ratio for a number of fluids in nature. 
T c denotes the critical temperature at the endpoint of the liquid-gas transition for water and helium, 
the supcrfluid transition temperature for ultracold Fermi gases and the deconfinement temperature for 
QCD. For the water and helium data (dotted lines) we refer the reader to [34]. The dashed curves are the 
expected theoretical curves for QCD (red) and ultracold Fermi gases (green), from [30, 32] and [35, 36] 
respectively. The solid red square denotes the upper bound ^ ~ 2.5 (^) for the QCD quark gluon plasma 
found in [37], while the open red squares denote the upper limits found in the lattice analysis of [65]. The 
dashed horizontal line is the universal ratio ? = t— ■ Similar plots can be found in [12, 38]. 

at LHC and RHIC, and may also help in finding the location of the critical point. However, 
at the moment most hydrodynamical simulations of the QGP assume that ^ is a constant, and 
therefore insensitive to temperature. The question of the possible relevance of a temperature- 
dependent 2 on the collective flow of hadrons in heavy ion collisions has been investigated in 
a number of studies [5, 39-41], which thus far have focused mostly on qualitative effects. The 
results of [39] seem to indicate that at LHC energies elliptic flow values are sensitive to the 
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temperature behavior of 2 in the QGP phase, but insensitive to it in the hadronic phase (with 
the results reversed at RHIC energies). 

Motivated by the potential sensitivity of elliptic flow measurements to thermal variations of 
2 at the energies probed by LHC, here we would like to initiate a systematic study of the flow 
of the shear viscosity as a function of temperature in the context of holography. In this paper 
we will restrict our attention to theories with a vanishing chemical potential, described simply 
by gravity coupled to a scalar field in the presence of higher derivative corrections. Although we 
will not consider string theory embeddings of these backgrounds, the higher derivative terms we 
examine are generically expected to show up as corrections to the two-derivative Einstein-scalar 
action viewed as an effective field theory, and in particular to most effective actions derived from 
string theory. Moreover, our result for 2 J n (18) is completely generic and applicable to any 
top-down model which contains a truncation to an Einstein-scalar theory. 

Thanks to the presence of a non-trivial scalar field profile, the higher derivative terms in our 
theory will generically generate a temperature flow for the shear viscosity, as already expected 
from [20, 21]. In particular, the temperature dependence of 2 will be dictated by the shape of the 
scalar potential and of the couplings of the scalar to the higher derivative interactions. Given an 
explicit expression for - in terms of the latter, it is then straightforward to discuss the existence 
of local minima for ^ as a function of temperature. As we will see, for potentials which confine 
quarks at zero temperature it is also possible to determine criteria for the existence of global 
minima. In particular, the requirements that - approaches its high-temperature value from 
below (dictated from asymptotic freedom) and that the zero-temperature theory is confining, 
are enough to fix the signs (and ranges) of the couplings of the higher derivative terms, and 
guarantees the existence of a global minimum for at or above T c . Thus, we have been able to 
give a holographic argument for the existence of a minimum for - in this class of models, along 
with a geometric interpretation for it, complementing existing field-theoretic 4 arguments [42]. 
Finally, although our analysis will be more general, we will focus mostly on 'phenomenological' 
models engineered to reproduce some of the qualitative and quantitative features of QCD. 



4 The authors of [42] tie a lower bound on the shear viscosity to the validity of second-order hydrodynamics. 
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The rest of the paper is organized as follows. Section II describes our setup and outlines how 
to obtain ^ for higher derivative corrections to dilatonic black brane solutions. Our main results 
are presented in Section III, where we discuss the temperature dependence of 2 i n theories with 
a non-trivial scalar field profile, first by considering a toy model consisting of an exponential 
scalar potential, then focusing on QCD-like phenomenological models. We also comment on 
qualitative features of ^ in theories that are confining, and identify generic criteria for the 
existence of a minimum for ^ a s a function of temperature in such setups. We summarize our 
results in Section IV, where we also comment on the challenges involved with determining more 
precisely the structure of -, and in particular how it flows with temperature. 



II. SETUP 

In order to generate a non-trivial temperature dependence for 2 at zero chemical potential 
we will consider backgrounds with a scalar field coupled to a higher derivative theory of gravity. 
The action we consider is of the form 

S = —V f d 5 *V=9 \R - 2(V$) 2 + V($) + t 2 P G($) R^paR^} , (1) 
lovrG 5 J L J 

under the assumption that the coupling f3 of the higher derivative terms is perturbatively small 
and G{Q) is an arbitrary regular function of For most of the paper we will focus on the 
exponential case, G(&) = e 7 *, however most of the analysis is straightforwardly generalized to 
arbitrary G{§). The dilaton potential V(<fr) is assumed to have either a minimum at some value 
$ = $o or a run-away behavior V{&) —> const as $ — > — oo (as in [43, 44]). Although in the 
above we have only considered the higher derivative correction R fJjL/pa R M ' upcT , this is in fact the 
only term which contributes to 2 at this order in the derivative expansion. The computation of 
2 in theories with higher derivatives has been studied in great detail, and here we review briefly 
only the relevant aspects 5 . 



For reviews of higher derivative corrections to - and further details of the computation we refer interested 
readers to [13-15] and references therein. 
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A. Extracting the Shear Viscosity to Entropy Ratio 

In the hydro-dynamic approximation to near-equilibrium dynamics, the transport coefficients 
of a finite temperature plasma can be extracted in a number of ways. The most straightforward 
method for computing the shear viscosity is based on the Kubo relation 

rj = - lim - Im (lj, k = 0) , (2) 

which reduces r\ to the low frequency and zero momentum limits of the stress tensor's retarded 
Green's function. Using the holographic dictionary, the relevant two-point correlation function of 
the shear stress tensor T xy can be read off from the effective action of a shear metric fluctuation 
h x v (t,u) ~ f d 4 k 4>k(u) e~ lut+,lkz , where we use u to denote the radial direction in the bulk. 
Expanding (1) to quadratic order in the modes <j) k , one finds the by now standard 6 form of the 
effective action for the shear fluctuation, 

+ D(u) 4> k <P-k + E(u) 4>i$'_ h + F(u) + S GH , (3) 

where the coefficients A(u),B(u), . . . ,F(u) encode information about the background solution, 
and Sgh is the generalized Gibbons-Hawking boundary term. 

After a number of subtle manipulations, the shear viscosity can be extracted from (3) and 
reduces [17] to the compact expression 

1 



V 



8vrG 5 



(4) 



U=U h 

evaluated at the horizon radius Uh, showing that it is given entirely in terms of horizon data 7 . 
Finally, the entropy density s is easily found by dividing by the (infinite) black brane volume 



6 See [45] for the original derivation, in the context of a' 3 corrections. 
The Kubo formula (2) can be shown [17, 46] to be equivalent to 



T) = hm - — , (5) 

u,^->0 IU <Pui,k=0(U) 

where Fife is the radial-momentum conjugate to <j>k- In the low frequency limit (as long as the boundary theory 
is spatially isotropic) the quantity inside the limit does not depend on the radial coordinate. As a result, it can 
be evaluated at an arbitrary value of u, and in particular at the horizon [17, 46]. 
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Wald's entropy formula, 

S = -2-ir / <fx\T-h— e^epa , (6) 

where h is the induced metric on the horizon cross section E, and the binormal to S. 

Prom our discussion above it is evident that ^ can be expressed entirely in terms of near- 
horizon data. Moreover, when working perturbatively in the coupling (3 of the higher derivative 
terms, ^ can be determined purely from the background solution of the two-derivative theory 8 . 
These two facts allow us to write ^ in terms of the parameters of a generic near-horizon non- 
extremal black brane expansion. Parametrizing the black-brane solution to the leading order 
two-derivative action by 

ds 2 = -a 2 (u)dt 2 + c 2 (u)du 2 + b 2 (u)dx 2 , $ = ip(u) , (7) 

with the choice b 2 (u) = we can write down its near- horizon expansion by assuming a first 
order zero in gu and a corresponding first order pole in g uu , 

a(u) 2 = a (l - u) + ai(l - uf + a 2 (l - u) 3 + ... , 

b(u) 2 = &o(l + (!-«) + ■■■), 

c(u) 2 = c (l - n)" 1 + ci + c 2 (l —u) + ... , 

ip(u) = ip h + (pi(l - u) + </> 2 (l - u) 2 + ... . (8) 
The shear viscosity to entropy density ratio is then of the simple form [20] 

l = (G(<p h ) + 2^i G'(cp h ) )) , (9) 

S 47T \ C / 

and is only sensitive to the parameters {cq, iph, yi} of the near-horizon expansions (8). This 
concludes the derivation of 2 for the theory described by (1). In the remainder of this section 
we will discuss the class of dilatonic black brane solutions we are interested in. 

8 Because of the universality of -, 0(/3) corrections to the background geometry lead to order 0(/3 2 ) corrections 
to 2. 
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B. Shear Viscosity of Dilatonic Brane Solutions 

As we mentioned above, the backreaction of the higher derivative terms on the background 
solution does not affect ^ to linear order in perturbative parameter j3, which is the order that we 
are interested in here. Thus, we are going to focus on black brane solutions to models described 
by the two-derivative action 



l — J d 5 x^ [r - 2(v$) 2 + , (io) 



16vrG 5 

and neglect curvature corrections. Although analytic black brane solutions are not known for 
generic choices of the potential, to extract - knowledge of the near-horizon behavior is enough, 
given our prescription (9). For this calculation we find it more convenient to introduce a new 
radial coordinate r and parametrize the black brane ansatz as 

ds 2 = f~ l {r) dr 2 + e 2A ^ (dx 2 - f(r) dt 2 ) , $ = $(r) . (11) 

From the expression (9) it is clear that we will need to calculate the near horizon values of the 
metric functions c(u) and ip(u) in (7), and the relationship between T and <ph, which will make 
the temperature dependence of (9) explicit. We will instead calculate the near horizon values 
of the functions 3> and / above and perform a change of variables at the end to express (9) in 
terms of physical quantities. 

For this study, it turns out to be particularly convenient to adopt the phase variables method 
developed in [47], which we briefly review here and in Appendix A. This is a quick and efficient 
way to obtain the thermodynamic properties of the diatonic branes. In place of solving the full 
5th order set of Einstein's equations, one only needs to solve two first order differential equations 
for the so-called "phase variables" that are defined from the metric functions as 

C $' If 

= Y ^^ih- (12) 

The constant Q depends on the normalization of the scalar kinetic term in (1) and dimension- 
ality. In our case it is fixed to be £ = \/8/3. Clearly, these functions are invariant under 
reparametrizations of the radial coordinate. Physically, one can interpret the boundary values 
of these variables as the thermodynamic energy and the enthalpy of the system. Furthermore, 



11 



their horizon expansion is completely determined in terms of the dilation potential V(3>/i) an d 
V'(&h) by the requirement of regularity. The behavior of the metric functions near the horizon 
is also determined in terms of these two quantities. The details of the calculation are explained 
in Appendix A and here we present only the final results. 

We parametrize the near-horizon expansion of the metric and scalar field of the black brane 
ansatz (11) at r = as follows, 

A(r) = A h + A 1 (r-r h ) + --- (13) 
f(r) = h(r-r h ) + --- (14) 
$(r) = $ fc + $ 1 (r-r h ) + -- - . (15) 

We can now make use of the phase variables method to obtain explicit expressions 9 for fi,A\ 
and $i in terms of the physical parameters in the system, 

fl = -M p (4vr)f^, (16) 

with C given by (A15). As can be seen from (9), the only near-horizon parameters that are 
needed for extracting - are cq, the leading order coefficient in the expansion of c(u), and ip^ and 
(pi, the first two terms in the scalar field expansion. By performing a change of radial coordinate 
we can relate the two near horizon expansions, and write {(fh, fl, Co} in terms of {A\, fi, $h, $1} 
as 

We now have all the ingredients needed to apply the near-horizon ^ prescription (9) to the 
generic black brane expansion (16)-(17) we just found. As expected, we find that the deviation 
from the universal result ^ = is controlled by the shape of the potential, the horizon value 



9 In particular, we use eqs. (A2), (A3), (A7), (A14) and (A16) in Appendix A. 
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of the scalar field as well as the coupling /5 of the higher derivative term, 




(18) 



and for the specific choice of G(<3?) = e 7 *, which we will focus on for the most of this paper, 




(19) 



We emphasize that this expression is completely general, and applies to any asymptotically AdS 
solution to (1). Moreover, we can already anticipate that ^ will generically be temperature 
dependent, thanks to the presence of a non-trivial scalar field profile, as already seen in [20, 
21]. To better understand the meaning of (19), we recall that the flow of - as a function of 
temperature is mirrored, in the bulk, by the change of the near-horizon geometry of the solution, 
as the horizon radius varies. For the class of holographic constructions we are considering here, 
it is the scalar field profile which is responsible for introducing an additional scale A in the theory 
(in addition to temperature), and thus breaking the conformal symmetry away from the UV. 
As a result, we should think of <&h, the horizon value of the dilaton, as tracking the temperature 
of the system 10 . Thus, (19) can be expressed entirely in terms of temperature by finding the 
precise relationship between <3?h and T. As usual, the latter can be determined from the metric 
by demanding regularity at the horizon. We review this is Appendix A. 

Finally, we would like to point out that in the special case of a non-dynamical scalar field 
and a constant potential, we recover the standard result for an AdS black brane in pure gravity 
with curvature corrections [48, 49], which is well-known to give rise to a constant correction to 
the universal 2 = result. More interestingly, there are special choices of non-trivial V(3>) for 
which the temperature dependence disappears, as we will see more explicitly below. 



In this section we examine the temperature dependence of 2 in various holographic setups. 
We start by looking at a toy model consisting of a simple exponential potential, and continue with 

10 More precisely, it will track the dependence on T/A, where A is the new scale in the system. 



III. RESULTS 
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more 'phenomeriological' constructions designed to mimic QCD and in particular the physics of 
the strongly coupled quark gluon plasma. We will conclude by remarking on generic, qualitative 
features of ^, including a discussion of the existence of minima as a function of temperature. 

A. A Warm-up Example: Chamblin-Reall Black Brane 

As a warm-up example, we will work out the Chamblin-Reall (CR) black-hole solution [50] 
by making use of the phase variables formalism. The CR brane is a solution to (10) with the 
single-exponential potential 

*W = ^e°* (20) 

where Vq is a positive dimensionless constant and t defines a length-scale in the background. We 
choose a > 0. With this convention 11 , in the zero temperature geometry the scalar approaches 
<]? — > — oo on the boundary and <!> — > +oo in the deep interior. In the black-brane geometry, 
which is what we are ultimately interested in, $ runs from — oo on the boundary to a constant 
value <£/j at the horizon. Thus, anywhere outside the horizon < <3>/j. 

Although the potential (20) does not admit an AdS minimum as the scalar field approaches 
the boundary, the CR brane - for a special value of a - can be obtained from a dimensional 
reduction of pure gravity plus a (negative) cosmological constant in six dimensions. In fact 
our action (1), with the potential choice (20), can be obtained via a U(l) reduction from the 
following six-dimensional mother theory 

£ = R + 2A 6 + pR^peRTf" , (21) 

for the special case 12 of 7 = —a = —y / 2/3 (see Appendix B for details of the reduction). This 
special parameter choice - for which the CR solution can be uplifted to a pure AdS black brane 
in six dimensions - will play an interesting role in the behavior of -, as we will see later in this 
section. For now, however, we will keep {a, 7} completely arbitrary. 



11 In the two-derivative theory, the sign of a can be altered by the transformation $ — > — $. 

12 More general dimensional reductions [51] may yield additional values of 7, a. 
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The CR solution corresponds to the fixed-point of the X-equation (A5) with 

ct 

X($)=x = -—, (22) 

where recall ( = i/8/3. Substituting this into the Y-equation (A6), the solution becomes 

Y($)= ,}~ X ° , K ={-a- x l). (23) 

From (A8) and (A9) one can then reconstruct the full background as a function of <J>, 

A(*) = A e +£-(*-* e ), /($) = 1 - <■«(**-*) , (24) 

where <3> c corresponds to a cutoff surface, as explained in the appendix. One can now use (A12) 
and (A13) to obtain the entropy and temperature as a function of <E>/j, 

T($ h ) = T e^^**, 5(* h ) = S e^*\ (25) 

where we defined 

To ss ^ e*-4» + *»»*- , S„ S S(0) = '"^^ . (26) 
The temperature dependence of the entropy can be read off from (25), 



S(T) oc r 1_4x o . (27) 

Combined with the first law of thermodynamics, this fixes the temperature scaling for the free 
energy of the system 

— 3 



Foc-T 1 " 4 ^. (28) 

In both equations (27) and (28) the proportionality constant is positive. Notice that the free 
energy F is always negative definite and never crosses zero - there is no phase transition in the 
system. Moreover, the requirement that the specific heat of the system C v = TdS/dT is positive 
definite (for thermal stability) constrains the value of xq, resulting in Xq < 1/4. Using (22) this 
can be translated into a condition on a, the exponent in the dilaton potential (20), 



0<a< x /-. (29) 
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Together with this condition, equations (27) and (28) essentially determine all of the thermo- 
dynamic properties of the system. The fact that - to work out the thermodynamics associated 
with this background - we only needed to solve two first order differential equations, (A5) and 
(A6), and not the full system of Einstein's equations, demonstrates explicitly the advantage of 
using the phase variables method. 

We are now ready to calculate the shear viscosity of the CR black-brane solution. Although 
we are particularly interested in the case of 7 = —a = —^2/3, for which the CR solution 
comes from a U(l) reduction of a pure AdS black brane in six dimensions, here we write down 
a 'formal' 13 expression for 2 for arbitrary values of a and 7. Using (20) and the corresponding 
expressions for the temperature and entropy (25), one finds 



T] 1 

S 4-7T 



— 6a(7+a) 

i-Vofi ^vn-™- 



(30) 



3' V 4 J \T 0j 

Recall from (29) that a 2 < 8/3, and therefore the sign of the power of the temperature in (30) 
depends on whether 7 < —a or not. In either case, 7 is a monotonic function of the temperature, 
and whether it decreases or increases relative to j- depends on the sign of /? as well as the range 
of 7. Furthermore, if we require the zero temperature limit to approach the universal 7- result, 
we must impose 7 < — a. 

Interestingly, the temperature dependence of the shear viscosity to entropy ratio disappears 
in the two special cases: 

— $ 

1. V oc e 3 T . In this case not only the T-dependence disappears, but also - resumes its 
universal value j-, despite the presence of higher derivative corrections. 



2. V oc e~ 7 *. For the special case of a = y/2/3, ^ takes exactly the same value as in the six- 
dimensional AdS Schwarzschild black hole. This fact can be understood by reducing the 
AcISq Schwarzschild solution on S 1 . It inherits the scale symmetry of the parent solution 
in six-dimensions, hence leading to absence of the T-dependence. See Appendix B for 
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Although our near-horizon - prescription was obtained specifically for geometries that are asymptotically AdS, 
it may be possible to generalize it to backgrounds that contain asymptotically conformally flat radial slices, by 
appropriately taking into account holographic renormalization in such backgrounds as in [52]. 
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details of this calculation. Similar statements can be made for reductions from (d+l + n)- 
dimensions on an re-torus as in [53] , which for d = 4 would yield a = —7 = ^ 3(^+3) • 

It is possible that the cancelation of the correction to ^ in the former case may also be under- 
stood in terms of a dimensional reduction of a parent theory, but this time without a cosmolog- 
ical constant 14 . At the moment we don't have a more complete understanding of this case (see 
however [54]). 



B. Improved Holographic QCD 

1. ihQCD Background 

Next, we turn to the phenomenological models discussed in [43, 44, 53] and focus in particular 
on the setup of [47, 55]. These are phenomenological constructions in the sense that the potential 
V(3>) is determined purely by field-theoretic requirements, and is designed to capture some of 
the features of QCD while remaining reasonably tractable. In fact, in these models the dilatonic 
scalar can be identified 15 with the running 't Hooft coupling, A = N c gyM ~ e*, and the scalar 
potential V(3>) is directly related to the /3-function of the system, giving a holographic definition 
of the latter in terms of the background geometry. Thanks to this identification one can extract 
the UV and IR asymptotics of the potential from the small A and large A expansions of /3(A). 

In the UV (small A), the input for the behavior of V(3>) comes from perturbative QCD, 
i.e. from the requirement of asymptotic freedom with a logarithmic running coupling. More 
generally, one assumes that there is a dimension-four operator TvF 2 in the spectrum, dual to 
the dilaton in the bulk. The fact that the operator Tr F 2 is marginal in the UV then translates 
into the statement that the UV geometry is asymptotically AdS$, with logarithmic corrections. 
The details of the UV expansion for these types of backgrounds can be found in [43]. 

On the other hand, in the IR (large A) the potential is fixed by demanding linear quark 
confinement. On the bulk side this is implemented by adding a probe string in the geometry - 



We thank Blaise Gouteraux for a very interesting discussion on this point. 
Up to a multiplicative factor which does not affect physical observables. 
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dual to a Wilson-loop - and obtaining the quark-anti quark potential from the asymptotics of 
the string embedding in the IR [44, 56]. One finds that linear confinement requires the dilaton 
potential to have an asymptotic expansion in the IR of the form 

V($) -> $ p + --- , (31) 



where we show only the leading large $ behavior 16 , with y/8/3 < Q < a/16/3 and P > 0. The 
parameter choice which fits the available zero and finite temperature data best turns out to be 
Q = a/8/3 and P = 1/2 (see [57]). This particular choice is also well motivated by the fact 
that it exhibits desirable qualitative features such as a linear glueball spectrum, screening of 
the magnetic quarks [44], and more recently the scaling behavior of the interaction measure in 
temperature [58] 17 . 

More specifically, defining A = exp(y / 3/2$), an example of a potential with the correct UV 
and IR asymptotics is of the form 



£ 2 V(\) 



4 



1 + A + ViAs 



log(l + V 2 \l + V 3 X 2 ] 



1/2 



(32) 



12 

where the value at A = sets the UV AdS scale I. The remaining parameters in the potential 
(32) can be fixed by matching the scheme-independent /3-function coefficients of large N QCD, 
the lowest glue ball mass and the latent heat at T c [57]. Although we will analyze explicitly 
the model with the potential given in (32), we emphasize that our qualitative results will only 
depend on the fact that V —> const in the UV and that it is a confining potential in the IR. 



2. Thermodynamics 

Before discussing the temperature dependence of the shear viscosity in this setup, we would 
like to summarize the basic thermodynamic properties of the system. These models exhibit a 
first order confinement-deconfinement transition at some critical temperature T = T C . Below T c , 
the dominant phase is the confined phase that corresponds to a thermal graviton gas background. 

16 The difference between this definition and the originai one given in [43] arises from different normahzations of 
the dilaton kinetic term. 

17 See [59] and references therein for a criticism of ihQCD models. 
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FIG. 2: Typical plots of the black-hole temperature (a) and free energy (b) as a function of the horizon 
position rh, in a confining background. The temperature features a minimum at r m i„, which separates 
the large black-hole from the small black-hole branches. The locus given by F(r c ) = 0, with r c < r m i n , 
corresponds to the phase transition point, T = T c . 

On the other hand, above T c one has the deconfined phase corresponding to a big black-hole 
background. We note that there also exists a third phase above a certain temperature T m j n , 
where T m j n < T c , which is sometimes referred to as a small black-hole. While the horizon of the 
latter is deep in the interior, the big black hole has its horizon closer to the boundary. 

The presence of the two types of black- holes is apparent from Figure 2(a), which shows T 
as a function of the horizon location. It is also clear from the figure that there are no black 
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FIG. 3: Black hole free energy as a function of temperature. The lines "B" and "S" denote the big and 
the small black-hole solutions, respectively. At T = T m i n the free energy of the two solutions is the same. 
The free energy of the thermal gas phase is set to F = 0. The direction along which the horizon location 
increases - which is represent by A^, the horizon value of the dilaton - is also shown in the figure. 

hole solutions for T < T m i n . Although the small black- hole is always sub-dominant and has a 
negative specific heat, its presence will turn out to be important to understand certain properties 
of - in the following discussion. Figure 2(b) shows the variation of the free-energy density F 
as a function of the horizon radius. Finally, by combining Figure 2(a) and Figure 2(b) one can 
parametrically solve for the free energy as a function of temperature. The result is sketched in 
Figure 3, which also summarizes the phase structure of the system as the temperature is varied. 

3. Shear Viscosity of ihQCD 

Given the general formula (19) and the dilaton potential (32), it is immediate to obtain 2 as 
a function of the scalar at the horizon. However, physically one is interested in having - as a 
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function of temperature, rather than Conversion from the latter to temperature is done by 
using either (A12) or (A16), after solving Einstein's equations numerically for the background 
functions. Figure 4 shows the results of this calculation, for the choice of potential (32) which 
gives the best fit to the available lattice data [57]. We note that the behavior of T as a function 

T 
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FIG. 4: Temperature as a function of A/j = eV^/^*' 1 in the ihQCD background. The labels A c and A m 
refer to the value of Xh at the critical and minimum temperature, respectively. 

of the horizon location is indeed of the form of Figure 2(a), with the minimum separating the 
big black hole (on the left) from the small black hole branch (on the right). 

Combining T(A/J with the analytic expression (19) and the potential (32), one can now plot 
2 as a function of T/T c , for a choice of the parameters (3 and 7. Unfortunately, we were unable 
to constrain the possible range of 7} with the available data, but instead chose representative 
values (recall however that we want /3 to be small enough so that the curvature corrections in 
(1) remain perturbative) . Depending on the choice of the couplings {/J, 7}, 2 w [\\ then display 
different qualitative behaviors. 

Two interesting fiducial cases - which are representative of the behavior of the viscosity 
in a large portion of the phase space - correspond to taking (3 > 0, < 7 < \/2/3 and 
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FIG. 5: Plots of - for theories with the scalar potential (32) and parameters chosen to give the best fit to 
lattice data [57]. In (a) we have taken ft = 0.1 and 7 = \/3/2, while in (b) ft = —0.01 and 7 = 6/5-^/3/2. 
Figure (a) shows that ^ has a small minimum near T w 1.8T C . We emphasize that the size of the variation 
in 2 can be reduced without changing its qualitative behavior, by taking ft sufficiently small. 

function of T which, for the parameter choices made there, appears around T ~ 1.8T C . In the 
high temperature limit ^ approaches the universal value from below. The behavior in the 
second case, depicted in Figure 5(b), is quite distinct. There ^ increases monotonically with 
temperature above T c , increasing indefinitely as T — > 00, just as in perturbative QCD. There is 
no local minimum appearing in the range T > T c , and therefore 2 acquires its minimum value at 
T = T c . In both cases, we know from field theory studies of the hadronic phase that ^ increases 
monotonically for T < T c , as one probes lower and lower temperatures. Thus, in Figure 5(b) we 
expect T = T c to correspond to the global minimum of the function ^(T). 

In these two examples, the qualitative behavior of the viscosity is different not only near 
T c , but also in the high-T regime. We should note, however, that our holographic results can 
only be trusted up to a certain T max , above which the perturbative higher derivative expansion 
breaks down 18 . In particular, we won't be able to trust an arbitrarily large deviation from the 



A potentially useful way of characterizing this break-down is by requiring i 2 fte 1 ^ R <C 1, and evaluating this at 
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universal ^— result. Other interesting qualitative behaviors are also possible, in the remaining 
range of 7}. For example, a mixture of the two aforementioned cases arises when /3 < and 
— y^8/3 < 7 < 0. In this case ^ exhibits a local minimum above T c just like in the first case 
above, but its high-T behavior is the same as that of the second case. We will provide a more 
detailed derivation of the different qualitative features of ^ as a function of temperature in the 
next section, where we will also discuss the boundaries of the various regions of phase space 
which lead to the distinct - behaviors. 

s 

C. Qualitative Features of the Shear Viscosity for Confining Backgrounds 

We would like to conclude this section by exploring some of the qualitative features of the 
temperature dependence of ^ from a more general point of view. For concreteness we will restrict 
our attention to dilaton potentials that exhibit confining IR asymptotics (as <1> — > 00) as in the 
ihQCD case of the previous section, 

^^e«¥|... (33) 

and AdS asymptotics in the UV. Furthermore, we will consider the following two possibilities 
for the behavior of the potential in the UV: 

a. V -> || + ve w * + • • • , as $ -00 

b. V -f^- ±m 2 $ 2 + • • • as $ -> . 

The first case (under the assumption that w > to ensure AdS boundary conditions) is precisely 
that of ihQCD-type backgrounds, where the dilaton is massless and corresponds to a marginal 
deformation by the dimension- four operator TrF 2 . The second case describes a massive dilaton 
of mass m, and corresponds to a deformation of the UV conformal theory by an operator of 

the horizon. This yields the relation |.£ 2 /3e 7 * ,r V(<l?ir) <C 1, from which one can compute T max , corresponding 
to the point at which the correction becomes 0(1). However, making this relation more precise is beyond the 
scope of this paper. 
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scale dimension 19 

Thus, relevant deformations correspond to m? < 0, and the Breitenlohner-Freedman (BF) bound 
is given by m?i 2 = — 16. 

In holographic constructions of the type we are considering, the flow of ^ as a function of 
temperature results from the way in which the near-horizon geometry changes as the horizon 
radius varies (it comes from sampling the phase space of the possible solutions to the theory). 
In our setup, the horizon value <&h of the scalar field tracks the temperature of the system. 
High temperatures then map to <&h — > — oo in the a) type potentials, and &h — > in the b) 
type potentials we have just discussed, whereas 'low' temperatures correspond to the region just 
above T ~ T c . We emphasize that the geometry which is relevant in this entire temperature 
range T > T c is that of the big black-hole discussed in Section IIIB 2. Thus, the high-T behavior 
of 2 will correspond to the large horizon limit of the big black- hole geometry. Below T c , the 
thermal gas phase dominates, and we won't attempt to compute ^ directly in that regime (we 
will discuss the subtleties involved with such a calculation towards the end of this section). 
Instead, we will extract information about the structure of ^ - and in particular any minima it 
has as a function of temperature - by exploiting the existence of both the big and small black- 
hole branches. Thus, even though the small black-hole is not of direct physical interest - it is 
not thermodynamically favored - it is still (indirectly) useful for probing -. 

We start by working out the viscosity to entropy ratio in the regime which corresponds to 
the small black-hole branch. In particular, we work in the &h >> 1 limit which describes an 
extremely small horizon radius (the far right of Figure 4). In this regime, we can extract ^ by 
plugging the IR potential (33) into our general formula (19), which gives 

^+"tK : rH e(,+e,t ** f+ '"' (35) 

19 Note that the normalization of the kinetic term for our scalar differs by a factor of four from that of more stan- 
dard AdS/CFT conventions, leading to a slightly different relation between the scalar mass and the conformal 
dimension of the dual operator. 
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where we are omitting terms that are subleading in the large <J? limit. Given that the small black- 
hole is not thermodynamically favored, we emphasize that the only crucial piece of information 
we need to extract from this expression is whether ^ increases or decreases as <&h grows larger 
and the small black-hole shrinks. We can see from (35) that as <3?f t grows, the correction to 2 
tends to become 20 large and positive for {/? > 0,7 > 4/3Q} and {j3 < 0,7 < 4/3Q}, while it 
tends into the large and negative direction in the opposite case (assuming in all these cases that 
7 + Q > 0). On the other hand, when 7 + Q < the sign of the exponential term changes, and 
the correction term vanishes in the extreme small black hole limit, from above (below) when the 
product [3 — 1^ is positive (negative). 

We note that in all of this analysis f3 is assumed to be perturbatively small 21 , \/3\ <C 1. Of 
course when the higher derivative correction (35) becomes 'too large,' the result can no longer 
be trusted - for consistency one would have to include not only the four-derivative term in 
the action (1), but all other higher derivative corrections as well. However, our expression for 
- above will still be useful for the following two reasons. First, in many phenomenological 
potentials of the type we are interested in, the leading IR behavior (33) doesn't necessarily set 
in at a very high value of &h- Thus, in these cases there will be a region of validity for (35), 
which can be made larger and larger by choosing j3 (or 7) smaller. Secondly, as it will become 
clear below, what we are really after is not the precise value of (35), but rather whether 2 tends 
upwards or downwards on the small black-hole branch. 

Next, we move on to discussing the high-temperature behavior of - realized on the big black- 
hole branch 22 , which maps to §>h — > —00 and <J>/j — > for type a) and type b) potentials, 
respectively. We analyze the two types of UV asymptotics for the potential separately. In case 
a) one finds, for negative and large &h (the black hole horizon is close to the runaway AdS 
region) 

2 -+ -L [1 - 8/3^ + • • -1 . (36) 
s 4ir 

20 Analysis of the sub-leading terms shows that it becomes large and positive even when 7 = 4/3Q. 

21 To ensure the absence of additional degrees of freedom associated with the higher derivative corrections in (1). 

22 We are now probing the near-horizon geometry of the big black-hole solution, in the large horizon limit. 
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Thus, in the high-T regime 2 attains its universal ^ value, for 7 > 0. Note that it approaches 
this value from above (below) for (3 < (f3 > 0). On the other hand, for 7 < 0, ^ tends upwards 
(downwards) for /3 < (/3 > 0). We note that in the 7 < case the result is only trustable up 
to a certain (negative) value of above which other higher order derivative corrections should 
also be taken into account. Thus, the expression for ^ is only reliable up to a certain T max . 
For the b) type potentials, i.e. with a true AdS minimum at <I> = 0, one finds instead 

^-3-(l-8/3)-^(l6 + m 2 £ 2 )<^ + ... , <& & -K). (37) 

First of all, we learn that 2 takes the value one expects 23 from a curvature-squared correction 
(with no dilatonic scalar coupling) in five dimensions [48]. Furthermore, recalling that the BF 
bound is m 2 l 2 > —16, we see that ^ approaches its constant high-T value from above (below) 
for /3j < (7?7 > 0;. 

There are a number of generic qualitative statements one can make about the running of ^ 
with temperature for these confining backgrounds: 

a. Divergence of : First, it is known [47] that for any potential that confines quarks 
at zero temperature, there exists a temperature T m j n (below T c ) in the finite temperature theory 
where the small and big black-hole branches meet, as shown in Figure 2a). This implies that 
diverges at this point and that 2 is double- valued above T m j n . As a result, although - itself 
is finite at T m j n , by the chain rule, its derivative gjf f = ^§r~M^^s w * n diverge there 24 . This is 
a generic feature of - that holds for any holographic background corresponding to a large N 
confining gauge theory at zero temperature. 

This behavior is exemplified in Figures 7 and 8, which show that the divergence occurs at the 
meeting point of the big and the small black-hole branches. We should note that it is the thermal 
gas, and not one of the black-hole solutions, that is the thermodynamically preferred solution 

23 This conclusion is independent of the sign of 7, because of the $/j — > limit. 

24 One may worry that a diverging might reflect a regime that can't be trusted. However, we can keep 
corrections to 2 arbitrarily small by appropriately tuning the perturbative coupling j3. More importantly, the 
divergence of the derivative of 2 a t T min is a result of the double- valued nature of $h(T) - the quantity 
itself remains finite and small. Thus, given that we are only interested in qualitative features of 2 - whether 
it increases or decreases near T m i n , and the sign of its slope there - the divergence of 2 is not expected to 
present a problem. 
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at this temperature and so the divergence in does not actually correspond to a physical 
property of the shear viscosity. Even though it cannot be considered a physical characteristic 
of ^ knowledge of the divergence of at T m , n will prove useful in analyzing the functional 
properties of 2. In particular, it will provide an holographic explanation for why ^ assumes its 
minimum value at the transition point T = T C for confining gauge theories, as we will see shortly. 

There is also a particular case in which the divergence occurs on the border of a thermo- 
dynamically dominant branch and, as such, does correspond to a physical property of 2. As 
described in [47] and [60], a second order confinement-deconfinement transition corresponds to 
the case when T m , n and T c coincide. In this situation there is no small black-hole branch, and 
the derivative of 2 generically diverges precisely at the transition point T = T c . More precisely, 
^ diverges at T = T c and therefore one generically expects that a i so diverges, as long as 
remains finite at that point. This latter condition should be analyzed separately. 

As a concrete example of this case, we take the generic backgrounds considered in [60], which 
exhibit a second (or higher) order Hawking-Page transition at a finite temperature T c at which 
&h — > oo. In order to find the behavior of ^ around the transition point, one should first make 
sure that the perturbative approximation in the higher derivative expansion is valid as $/i — > oo. 
The scalar potential in the theories of [60] has the asymptotic behavior shown in (33) with P = 
and Q = a/8/3, together with exponentially decaying sub- leading terms as <3? — > oo, 

-> eV 7 !* (l + Ce~ K *) , k>0. (38) 

Here C > is a positive constant which does not play any role in the following discussion, and 
the value of k determines the order of the continuous transition. 

From the 2 expression (35) we see that we can only make reliable statements about the 
behavior of the shear viscosity as T — > T c (&h — > oo) in the particular case 7 < — Q = —a/8/3 
(otherwise our perturbative approximation is guaranteed to break down). For this choice of 
parameters one then recovers the universal result ^ -> ^ as T T c . It is straightforward 
to work out the behavior of the derivative of ^ by combining 2 and The first can 
be calculated straightforwardly from (35). The latter follows from the sub- leading analysis 
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performed in [60], which gives the following relation between temperature and 

r($ A )^T c + Ce-*, *h->oo, (39) 

where C > 0. Finally, combining these ingredients one obtains, 

A WW§ * fc -oo, (40) 
ai s 

where the sign of the proportionality constant is positive (negative) for f3 < (/3 > 0). Thus 
we arrive at the conclusion that, although r?/s remains finite and reaches its universal ^— value 
at the continuous transition, its derivative diverges for —a/8/3 > 7 > — ^8/3 — k, and it 
vanishes for 7 < — a/8/3 — k. We emphasize again that our perturbative approximation breaks 
down for 7 > — -y/8/3, in which case one cannot reliably determine the behavior of ^ at T = T c . 
This concludes our aside on the special case in which the divergence occurs on the border of a 
thermodynamically dominant branch. In the rest of the paper we will assume that T c 7^ T m j n . 

b. Conditions for local minima: Our general expression (19) makes the discussion of 
the presence or absence of a minimum for - straightforward. Local minima are defined by the 
following two conditions 25 , 

d 77 d 2 11 , , 

-^-- = 0, ^2->0, 41 
d$ h s d<b z h s 

which, when re-instating the generic higher derivative coupling G(&), can be translated into 
analyticity conditions on G(3>) and the dilaton potential V(<&) using (19). In particular, extrema 
of - are determined by the locii &h at which 

jG'(<f>)V'($ h )-G($)Vm) m =0, (42) 

an expression which is valid for arbitrary potentials (which allow for AdS asymptotics) and 
scalar couplings G, and not just the particular setup of this section. In terms of G($) = e 7 *, 
the condition becomes 

^yV"($ h ) + (^ 7 2 - W($ h ) - 7 V(* h ) = . (43) 



25 Without loss of generality we assume that dT/d&h is finite at these minima. 
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Thus, in order for 2 to have an extremum at a certain temperature, a necessary condition 26 is 
that the dilaton potential V(<&) possesses at least one solution to equation (42). Clearly, for 
such a solution to correspond to a minimum one should further ensure that the second derivative 
is positive there. These two conditions provide a simple criterion for the existence of extrema, 
given an arbitrary scalar field profile in the class of theories (1). 

c. Presence of a global minimum: The discussion in the previous paragraph concerns 
general conditions for the existence of local minima. However, we can also ask what are the 
conditions for the existence of a global minimum, given a confining potential with the IR and 
UV asymptotics described above. To answer this question, it turns out to be useful to recall the 
behavior of 2 in the two opposite limits of an extremely small and extremely large black-hole. 
We will restrict our attention to the case in which the function -(T) is monotonic on the small 
black-hole branch, an assumption which is satisfied in all of the holographic backgrounds we 
consider in this paper. Our analysis will also make use of the fact that the derivative of 2 
approaches ±oo as T — > T m i n along the two black-hole branches, as we discussed above. We will 
divide the discussion into two cases: 

1. Discontinuous case: Here we discuss the situation in which 2 displays a global minimum 
appearing exactly at T = T c . If - is a monotonically increasing function of temperature 
on the entire big black-hole branch, reaching its maximum value as T — > oo from below, 
then clearly it will reach a minimum value (on the big black- hole branch) at T = T c , 
where the big black hole solution ceases to be thermodynamically favored. To work out 
the behavior of - in the range T < T c one would then need to switch to the thermal gas 
solution. Calculating 2 directly in this regime is technically very challenging. However, 
as outlined in the introduction, there are qualitative field theory arguments which tell 
us that - should start increasing again below T c , as the temperature continues to drop. 
Thus, we conclude that in this case T = T c corresponds to a global minimum for 2 This 
behavior is shown schematically in Figure 6, from which it is also clear that the derivative 
of - is discontinuous at T r . 



This is also sufficient provided that dT/d&h is finite at this point. 
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FIG. 6: A cartoon of the behavior of ^ in the big black-hole phase (to the right of T c ) and in the thermal 
gas phase (for T < T c ) for the case where the big black-hole phase is monotonically increasing above T c . 
We denote by T max the temperature above which the correction to j is no longer perturbative, and our 
approximation breaks down. 

Graphically it is clear that, in order for this to happen, on the big black- hole branch 2 
must approach its high-T value from below, and furthermore its derivative on this branch 
must approach — > +oo at T m j n , where the big and the small black- hole branches 
meet. Given our assumption that 7 is a monotonic function of T on the small black- hole 
branch, in order to satisfy the second criterion it suffices that ^ decreases as $^ gets larger 
(in the extreme $^ 3> 1 limit) on the small black-hole branch. By examining (35) and 
(36) we see that this requires 27 either of the following two cases: 

[i] > 0, < 7 < ^ , (44) 
[«.] P < 0, 7 < -Q . (45) 

27 Note that we need the derivative of 2 to be positive definite in the asymptotic high-T region (on the big 
black-hole branch). Since dT/d&h is negative definite on the big black-hole branch, we see from (36) that this 
requires /J7 > 0. 
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FIG. 7: (a) A cartoon of the case when (44) is satisfied, for the type a) potential, (b) A cartoon of the 
case when (45) is satisfied, for the type a) and type b) potentials labeled respectively by "case a" and 
"case b." In both figures T m denotes the temperature where the small (S) and the big (B) black-hole 
branches meet and T c denotes the temperature of the confincment-dcconfinement transition. The high 
temperature value of the viscosity to entropy ratio is denoted by 
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We plot these two cases schematically in Figure 7(a) and 7(b). Figure 7(a) shows the 
behavior of ^ on the big (B) and small (S) black-hole branches for a potential whose UV 
asymptotics are described by case a). In Figure 7(b) on the other hand we include both 
cases a) and b). 

2. Continuous case: Here the 2 curve has a global minimum on the big black-hole branch 
at some temperature T„ satisfying T„ > T c . Note that at that point, - and all of its 
derivatives are continuous. For this to happen, clearly - must increase as the temperature 
is raised above T„, approaching its high-T value from below. More importantly, 2 must also 
increase below T n , as the temperature gets closer and closer to T m j n . The latter condition 
translates into the requirement that the derivative on the big black hole branch obeys 
— > —oo at T m i n . Finally, it suffices that 2 tends upwards on the small black- hole 
branch 28 , in the &h 3> 1 limit. Investigation of (35) and (36) reveals that one needs 

[*■] P > 0, 7 > ^ , (46) 
[it.] P < 0, -Q < 7 < 0. (47) 

We provide schematic drawings of these two cases in Figure 8(a) and 8(b). Note also 
that (46) corresponds to the case plotted in Figure 5(a). Of course one should also incor- 
porate the fact that, in confining theories, the thermodynamically preferred background 
changes below T c . Thus, in order to be able to observe a global minimum located at some 
temperature T„, one should ensure that T v > T c , as shown in Figure 8(a). 

The four cases depicted in Figures 7(a), 7(b), 8(a) and 8(b) cover most of the physically interest- 
ing situations. Qualitatively distinct behavior (for example a maximum rather than a minimum) 
for - arises in the other ranges of (3 and 7, but we will not discuss it here. Finally, we should 
note that allowing for a more general scalar coupling G(§) to the higher derivative term [i.e. 
of the racetrack-type) would modify - and potentially significantly complicate - the analysis of 
the existence of a global minimum. We will not attempt this here. 

28 Recall that we are assuming that ^ is monotonic there. 
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FIG. 8: (a) A cartoon of the case when (46) is satisfied, (b) A cartoon of the case when (47) is satisfied. 

denotes the high-T value of -, T m denotes the minimum temperature where the small (S) and 
the big (B) black-hole branches meet and T c denotes the temperature of the confinement-deconfinemcnt 
transition. In Figure (b) we plot both the type a) and type b) potentials labelled by "case a" and "case 
b" respectively 
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1. The Thermal Gas Phase 

We would like to close this section by discussing briefly the main difficulties involved with 
calculating ^ G n the thermal gas background which describes temperatures below the decon- 
finement transition. We emphasize again that in our analysis we did not explicitly compute - 
below T c , but rather used knowledge of its behavior in the hadronic phase. The main difficulties 
involved with a direct computation for T < T c stem mostly from the need to consider 1/N 2 
corrections. 

To see why this is so, note that the thermal gas solution has no horizon, hence the 0(N 2 ) 
contribution to the entropy that would come from evaluating the action on the classical saddle 29 
vanishes. However, the finite temperature in the background generates thermal fluctuations 
of the graviton gas. Thus, the entropy in this case should be calculated by computing the 
determinant of fluctuations around the classical thermal gas saddle, and as such it is 1/N 2 
suppressed with respect to the 0(N 2 ) black- hole contribution. As a result, the entropy of the 
thermal gas is st g = 0(1). 

Now, let's move on to the calculation of r\ on the thermal gas background, which can be done 
by considering fluctuations of the transverse traceless gravitons - in particular, ^ can be related 
to the flux of transverse gravitons through the horizon, which is proportional to e 3Ah N 2 . On 
the other hand, the small black-hole background asymptotes to the (Lorentzian) thermal gas in 
the limit of vanishing horizon radius. This tells us that the leading 0{N 2 ) contribution to rj 
also vanishes, in the same way in which the 0(N 2 ) contribution to the entropy of the thermal 
gas vanished. In order to find the finite contribution one has to consider Witten diagrams with 
quantum loops, which lead to rjt g oc 0(1). Thus, although ^ in the thermal gas is of the same 
order as that for any black-hole background, in order to calculate it one needs to compute 
the 1/iV 2 contributions to both rj and s. Such a calculation would involve computing one-loop 
diagrams in supergravity. This would naively require knowing the spectrum of all supergravity 
fields which could run in the loops 30 and as such, knowing a genuine string theory embedding 



Note the identification 1/Gs oc N 2 in front of the action. Thus, whenever the on-shell action on a classical 

saddle is non- vanishing, automatically all of the thermal functions are proportional to iV 2 . 

Note, however, that one-loop contributions to hydrodynamic long-time tails do not require a detailed un- 
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of any particular model. This one-loop calculation appears to be a difficult task and we don't 
attempt it here. 

IV. DISCUSSION 

Before summarizing our results for the temperature dependence of the shear viscosity 
in holographic models, we would first like to discuss several properties of rj/s in the QCD 
plasma. Although the viscosity to entropy ratio in the QGP phase is expected to be very small 
and roughly comparable to ^— , its precise value is still largely uncertain. The collective flow 
observed at RHIC and LHC is analyzed by performing a Fourier decomposition of the particles' 
angular distribution, whose Fourier coefficients v n are sensitive to transverse momentum px- 
Restrictions on the possible range of ^ are then extracted from the data either by comparing 
the momentum-dependent elliptic flow coefficient V2{pt) with results obtained from viscous 
hydrodynamical calculations 31 , or by fitting the centrality dependence of the average pt 
integrated elliptic flow. Both approaches yield bounds on - which are somewhat close to the 
universal ^— value, with the analysis of [1] giving Air 2 < 2.5. 

Uncertainties of 2 

s 

A number of challenges, both theoretical and experimental, have to be overcome in order 
to determine ^ more precisely. Uncertainties on the experimental side come from non-flow 
correlations introduced, for example, by jets and resonance decays. On the theoretical side, 
uncertainties are due to various assumptions on the choice of initial conditions, initial state 
pressure gradients and event-by-event fluctuations. As an example, it was realized only recently 
that the assumption that higher order harmonics were negligible was a poor one - the shape of 
the fireball can fluctuate from one event to the next, even at a fixed impact parameter. The 
resulting irregular pressure gradients are not symmetric with respect to the reaction plane, and 

derstanding of the spectrum of supergravity fields [61]. Perhaps a similar simplification occurs for one- loop 
contributions to r\j s. 

31 More precisely, one needs a 'hybrid' code which combines viscous hydrodynamics of the QGP phase with a 
realistic model of the late hadronic stage (see e.g. [62]). 
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can in fact induce higher harmonic flow patterns. The elliptic flow direction and magnitude can 
also fluctuate event-by-event. Interestingly, the measurement of higher harmonics can play an 
important role in reducing the degeneracy between the shear viscosity and initial conditions, 
and appears to give tighter limits on ^. In fact, a combined analysis [63] of elliptic and 
triangular flow coefficients (v2 and 1)3, respectively) in Pb + Pb collisions at LHC was shown 
to favor a small shear viscosity ^ ~ 0.08, disfavoring the considerably larger value of ^ ~ 0.2 
(we refer the reader to [7] for a more detailed discussion of these issues). Although none of 
the models currently used describes perfectly all experimental flow data, the expectation is 
that - with the new high precision data coming from LHC - it will be possible to cut down on 
such uncertainties and arrive at a precision measurement of the transport properties of the QGP. 

The temperature dependence of ^ 

Another important issue for disentangling the physics of the QGP is that of the possible temper- 
ature dependence of -, which is the main focus of our analysis. Although most hydro dynamical 
simulations of the QGP assume a constant value for ^, a number of studies - mostly qualitative 
in nature - have begun to examine the possible relevance of temperature in LHC and RHIC 
heavy ion collisions [5, 39-41]. First results for the shear and bulk viscosity obtained via lattice 
QCD simulations were put forth in [64, 65], where it was shown that ^ remains close to its uni- 
versal ^ value for a range of temperatures not far above T c . On the other hand, kinetic theory 
results for QCD at weak coupling [32, 33] imply 2 > 1, for the QCD running coupling constant 
in the range a s < 0.25. Elliptic flow values at LHC energies may be sensitive to the temperature 
behavior of ^ in the QGP phase, although insensitive to it in the hadronic phase [39]. Moreover, 
it was argued in [28] that - should have a minimum at (or near) the QCD phase transition. 
This is expected because - increases with decreasing temperature in the hadronic phase [29, 30], 
while asymptotic freedom dictates that it increases with temperature in the deconfined phase 
- thus leading to at least one minima somewhere in the intermediate region. An interesting 
study of shear and bulk viscosities for a finite temperature, pure gluon plasma was performed in 
[41], where the authors used a phenomenological, quasi-particle model based on an effective ki- 
netic theory description. At large temperatures, the results of [41] reproduce parametrically the 
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transport behavior observed in perturbative QCD calculations. As the non-perturbative regime 
is approached, their analysis is in agreement with lattice QCD results, with a decrease of ^ with 
temperature and a minimum at T > T c . Independently of the assumptions behind particular 
choices of model and simulation schemes, it is clear that a more systematic understanding of 
the temperature dependence of ^ is an important ingredient towards a better descriptions of the 
dynamics of the QGP. 



The aim of this paper was to initiate a systematic study of the temperature dependence of the 
shear viscosity in the strongly coupled QGP (and gauge theory plasma more generally), in the 
framework of the holographic gauge /gravity duality. We were largely motivated by the potential 
sensitivity of elliptic flow measurements to thermal variations of ^, but were also interested in 
gaining further insight into the geometric interpretation of a non-trivial temperature flow for 
the shear viscosity. We restricted our attention to theories of gravity coupled to a scalar field in 
the presence of higher derivative corrections, 



under the assumption that the coupling /3 is perturbatively small 32 . Higher order curvature 
corrections are well-known to push 2 away from its universal ^ value. Moreover, once they are 
coupled to a non-trivial scalar field profile, they lead generically to a temperature flow for ^, as 
already seen in [21] and emphasized in [20]. 

The correction to 2 - which is parametrized entirely in terms of horizon data - can be 
expressed in terms of the potential and the horizon value of the scalar field, and for (48) 
takes the relatively simple form (see Section II for details of the derivation) 



A. Summary of Results 



C = R- 2(V$) 2 + V{§) + PG(§) R^ upa R^ pa 



(48) 



= h 1 + \ P [- G (®h)V{$ h ) + ^G'($ h )V\* h ) 



) 



(49) 



From now on we will set the AdS length scale I equal to one. 
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Note that thanks to this expression, determining the presence of local minima for 2 i s entirely 
straightforward. In particular, extrema are found by minimizing (49) with respect to 33 <3?h> and 
therefore correspond to the zeros of the relation (42). We emphasize that such a condition is 
rather generic, given the broad assumptions behind (49) - essentially the requirement of AdS 
asymptotics. 

Armed with the general expression (49), we have analyzed a number of holographic setups. 
For the majority of our analysis, we have restricted our attention to the special case of G($) = 
e 7 *, for which the viscosity to entropy ratio reduces to 



V 1 

s Air 



l + lpf-V($ h ) + l'yV'(* h ) 



(50) 



In this class of theories, the (horizon) quantity tracks the temperature of the system. As a 
result, the behavior of ^ as a function of temperature is controlled by the scalar field coupling to 
the curvature corrections - parametrized here by 7 and /3 - as well as by the specific functional 
form of the potential. 

As a simple consistency check of our result, we note that for the case of a non-dynamical 
scalar our expression (50) reduces to a constant and reproduces the standard result for an AdS 
black-brane in pure gravity with curvature-squared terms. On the other hand, an analytically 
tractable example which gives rise to non-trivial temperature dependence is that of a single 
exponential potential, V tx e a *. Black-brane solutions in this system have a temperature of 
the form T oc e 3o A 8 / k . With this choice of potential, our expression (50) implies a 
monotonic temperature flow for ^, whose precise structure is shown in (30) and is dictated by 
the range and signs of j3 and 7. Interestingly, when 7 = —a the viscosity to entropy ratio 
takes on a constant value again, with no temperature dependence. In particular, for the specific 
choice 7 = —a = a/2/3 the model can be obtained via a U(l) dimensional reduction of a six- 
dimensional theory of pure gravity with a negative cosmological constant, in the presence of 



R corrections. As expected, our correction to 2 in this case reduces to that for a pure AdS t 



6 



black brane, providing a rather non-trivial check of our result (50). The absence of temperature 
dependence in this case is then a simple consequence of the fact that the five-dimensional theory 



Assuming that dT/d<&h is finite. 
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inherits the scale invariance of the parent six-dimensional one (see also the discussion in [54] ) . 

The behavior of ^ becomes more interesting - and more relevant to the physics of the QGP - 
for theories which undergo a confinement-deconfinement transition at some critical temperature 
T c . In our analysis we have worked mostly with potentials constructed 'phenomenologically' 
by requiring their IR and UV asymptotics to match with expectations from QCD, as in the 
case of the Improved Holographic QCD (ihQCD) model discussed in Section IIIB. However, 
we emphasize that a number of our results hold more broadly, and can be easily generalized. 
Plots of the temperature dependence of - in the ihQCD model for fiducial values of {/3, 7} 
are shown in Figure 5. On the left, we note the presence of a small minimum slightly above 
T c , above which 2 approaches -h-. The behavior for T < T c was not computed directly, but 
rather extrapolated using the fact that ^ is expected to increase in the hadronic phase, as the 
temperature is lowered. The figure on the right, on the other hand, exhibits a minimum exactly 
at T = T c , with 2 growing indefinitely as the temperature is raised. A cartoon of the same type 
of behavior - expected to describe more general settings - is shown in Figure 6. 

For theories which confine quarks at zero temperature, the thermodynamic structure of the 
corresponding finite-temperature theory plays an interesting role in determining the behavior of 
2, and enables one to make generic qualitative statements about its temperature flow. In this 
class of models, the geometry which is favored thermodynamically for T > T c is that of a big 
black hole, while below T c the dominant phase corresponds to a thermal gas. However, a third 
solution appears for T > T m i n (with T m i n < T c ), which describes a small black- hole solution. 
For the ihQCD model, the existence of the two black-hole branches can be seen in Figure 4, 
with the minimum separating the big black hole (on the left) from the small one (on the right). 
Although the small black-hole is never thermodynamically favored - and is therefore not of any 
direct physical relevance - it can still be useful for probing - in the physically relevant regime 
T > T c . In fact, the coexistence of the small and big black- hole solutions for T > T m i n implies 
that - is double-valued there, and that its derivative diverges at the point T = T m i n where 
the two branches meet. These features can be seen for example in Figure 7(a). By examining 
the rough dependence of ^ on temperature on the two branches and the sign of as T m i n is 
approached, one can then reconstruct whether ^ exhibits a minimum in the T > T c range. 
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Following this strategy, in Section IIIC we have derived a set of geometric conditions for the 
existence of a global minimum for ^ in this class of models 34 . We have chosen to work with dila- 
tonic potentials which have confining IR asymptotics and exhibit two distinct behaviors in the 
UV, both consistent with AdS asymptotics - the type a) and type b) potentials corresponding, 
respectively, to a massless and massive scalar. The requirement - from asymptotic freedom - 
that ^ approaches its high-T value from below fixes the signs and ranges of the higher derivative 
couplings f3 and 7. Combined with the requirement that the zero-temperature theory is confin- 
ing, it guarantees the existence of a global minimum for - at or above the critical temperature, 
in a certain portion of the phase space. Furthermore, whether the (global) minimum is located 
exactly at T c or above it can be determined by looking at the sign of the divergence of at 
T m i n , the point where the big and small black hole branches meet. 

In Section IIIC we have provided a geometric classification of these two cases, which we 
emphasize are physically distinct. The parameter space of the ' discontinuous 1 case, describing a 
global minimum exactly at T c , is summarized in (44) and (45), and shown schematically in Figure 
7. The ' continuous 1 case, with a global minimum at some temperature T > T c , is described 
instead in (46) and (47), and shown in Figure 8. We note that both cases have (3j > (for the 
precise range of 7 we refer the reader to Section IIIC). The two types of UV asymptotics, cases 
a) and b), are also shown in the figures. The high-temperature behavior of ^, and in particular 
whether it approaches a constant value or increases indefinitely with temperature, is determined 
by the values and ranges of the /?, 7 parameters. Another generic result of our analysis is that 
when the confinement-deconfinement phase transition is second (or higher) order the derivative 
diverges or vanishes at T c , depending on the value of the parameter 7, while n/s remains 
finite and approaches its universal j- value 35 . 

One should note, however, that the holographic demonstration of the existence of a global 
minimum in confining backgrounds is incomplete. In particular, it is desirable to find an indepen- 
dent holographic reason for the requirement that /3j should be positive, needed for the presence 
of the global minimum in the continuous/discontinous cases discussed above. In fact, although 

34 Under the (mild) assumption that 2 ; s monotonic on the small black-hole branch. 

35 This result is valid as long as the perturbative approximation can be trusted. See section IIIC for details. 
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we always insisted on keeping the correction to ^ perturbatively small (by appropriately turning 
down /3), we relied on the fact that diverg es near T m in- Thus, one may worry that any 
conclusion drawn from using ^| in a regime where it is very large might not be valid. While 
this is a fair concern, we should note that it is which diverges, with ^5-^ remaining finite 
(and small). For completeness, it would also be desirable to perform the analogous ^ calculation 
in the hadronic phase (T < T c ), and derive holographically the monotonically decreasing behav- 
ior of - expected from field theory. Finally, we should caution the reader that a more generic 
higher derivative scalar coupling G(&) would give rise to a more complicated structure for 2 
(i.e. imagine allowing for racetrack-type terms) and might invalidate some of our arguments for 
the presence of a global minimum. Clearly, our analysis would have to be re-evaluated in such 
situations. Regardless, we would like to emphasize that - for a certain class of potentials, and 
for simple higher derivative couplings of the form G oc e 7 * - we have translated the issue of the 
existence of a minimum for ^ into a geometric one, and identified the holographic conditions 
which would guarantee its presence, complementing the field-theoretic arguments of [42]. 

In order for our analysis to be of more direct relevance to studies of the QGP, it is crucial 
to find ways to place restrictions on the allowed ranges of the couplings /3 and 7, whether 
with theory or experiment. This would cut down on the large phase space for the temperature 
behavior of and eliminate some of the model dependence inherent in holographic setups of 
this type - in particular when higher derivative terms are present. At the moment, however, 
this is very challenging. In the future, it might be feasible to constrain models by examining the 
effects of higher derivative corrections on the remaining transport coefficients (including those of 
second order hydrodynamics), and then simultaneously fitting to the available data. While this 
should be possible in principle, at the moment it presents a serious challenge - the bulk viscosity 
flows with temperature already without the need for higher derivative terms, and very little is 
currently known from holography about second order transport coefficients in the presence of 
higher derivative corrections. On the other hand, it may be possible to restrict the values of the 
parameters with more phenomenological considerations. As an example, at least in principle we 
can fit the UV behavior in the type a) potentials we discussed to that predicted by perturbative 
QCD, which would fix the values of both couplings f3 and 7. However, we should keep in mind 
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that our analysis is not valid to arbitrarily high temperatures (it breaks down when the higher 
derivative interactions are no longer perturbative), and therefore any direct comparison with 
UV physics should be taken with a grain of salt. We leave these questions open for future work. 
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Appendix A: The Method of Phase Variables 

For the study of dilatonic black hole solutions, it turns out to be particularly convenient to 
adopt the phase variables method developed in [47], which we briefly review here. With the 
black brane ansatz 



ds 2 = f-\ r ) dr 2 + e 2Mr) (dx 2 - f(r) dt 2 ) , $ = , 



(Al) 
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the Einstein and dilaton equations of motion that follow from the action (10) can be put in the 
following simple form of five first-order equations, 

d A_ = * -CtfJW* (A2) 
dr t 

d ± = * jr(*) e -</o**(*>* (A3) 
dr £Q 

\*L = -^(^e-Cr^, (A4) 



and 



where ( = and we have introduced the phase variables 

C If 

y($) "4/^' (A7) 

which are invariant under radial coordinate transformations [47] . Note that the metric functions 

A and / can be expressed in terms of X and Y by integrating directly ( A7) , 

C /"* d§ 

A(*) = A($ c ) + ^ , (A8) 

/($) = exp(cj*^d$) , (A9) 

where <£ c denotes a cut-off surface that plays the role of the regularized UV boundary 36 . 

The advantage of expressing Einstein's equations in terms of the phase variables {X, Y} 
should now be apparent - the system has been reduced to a coupled set of first order equations , 



36 At the boundary we require / — > 1, which fixes its dependence on <E> C . Moreover, as shown in [47], the physical 
observables in this system only depend on X and Y and are independent of <£> c . 

37 In order to solve (A5) and (A6) one has to specify a single boundary condition for each of the equations. After 
demanding a regular horizon of the form 

/($) = const, x ($h — $) near $ « $h, (A10) 

inspection of the equations (A9) and (A6) shows that the remaining integration constant is completely fixed, 
and near $ ~ $h one must require 

= "I* W + - $) ' Y($) = + ■ (AU) 
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(A5) and (A6), which greatly simplifies the task of finding analytic solutions. Moreover, all 
thermodynamic properties of the system are completely determined by knowledge of the two 
degrees of freedom X and Y as a function of <!>. In fact, in terms of the horizon value &h of the 
dilaton, the temperature and the entropy of the black brane are given by 

r(*h) = 4~ eA{ " h) ^-^ hx ^ d * , (A12) 

■e 3A ^. (A13) 



4Gat 

Combining these two expressions we find a simple relation between the phase variable X and 
the scalar potential, 

e -<^X W » =C ^ V ^ h)t (A14) 

where the proportionality constant C is given by 

and we have traded Newton's constant for the Planck mass, = (167TG7V)- 1 . Finally, the 
temperature can be directly related to the horizon expansion of the blackness function as follows, 

4vrT = - f (r h ) e A{Th) . (A16) 

These relations greatly facilitate the study of black brane solutions to the model described by 
(1), as we show in the main text. 

Appendix B: Chamblin-Reall from U(l) Reduction 

In this appendix we expand upon the particular case of the Chamblin-Reall black brane. The 
special case 7 = —a, for which 2 shows no temperature dependence, can be realized as a solution 
of a dimensionally reduced theory of pure gravity plus a negative cosmological constant in six 
dimensions. As such this solution can be uplifted to a pure AdS black brane in six-dimensions, 
thus explaining the temperature independent form for -. 
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To see this in more detail, we start from the following Lagrangian density in d+ 1 space-time 
dimensions, 

C = R + 2A d+1 + pR^STi" , (Bl) 

and perform a circle reduction by decomposing the metric as 

ds 2 d+l = e 2 ^ds 2 d + e 2<5< ^ 2 , (B2) 

where <p is a scalar 38 . The dimensionally reduced theory is then given by 

C d = R- a 2 (d-l)(d- 2)(9<A) 2 + 2A d+1 e 2 ^ + Pe~ 2 ^ \R a ^ s R aMS + ...] , (B3) 

where we fixed 5 = —{d — 2) a to ensure a canonical Einstein term, and omitted higher derivative 
terms involving R^ u and <fi, since they do not affect 2 Notice that, as expected, we have 
generated a dilatonic coupling to the Riemann-squared term. However - unlike in (1), where 
the parameter 7 is free - here the dilatonic coupling is completely fixed in terms of a. Thus, we 
have recovered the model of (1) with the CR potential (20), but only for the special parameter 
choice 7 = — q. 

We can now compute 2 for the dimensionally reduced theory. Restricting to d = 5 and 
rescaling the dilaton so that its kinetic term agrees with (1), the lagrangian becomes 

£5 = R-2(d(P) 2 + 2A 6 eV^^ + f3e~V^4>\ R ^/S Ra ^ s + _] ( B4 ) 

Comparing this to our action (1), we see that we have generated a scalar coupling ~ 
e 7<1 ' R f iupaR tJ,l ' pa with 7 = —y^2/3. The black-brane solution to the two-derivative theory is 
given in [66], 

ds 2 = e 2A (-f(r)dt 2 + dx 2 ) + (B5) 

f{ r ) 

e A = r 4 , (B6) 
f(r) = l-(y) 15 , (B7) 

(B8) 



f2/U - Aer " 



90 



38 Note that we have consistently fixed the Kaluza-Klein gauge field to vanish. 
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and has the following near horizon expansion 



? h (r-r h ) + ..., 

lnr h + 7^( r ~ r h) + -J 

4> h -^(r-r h ) + .... 



A{r) 



cj){r) 



(B9) 




This is precisely the result expected for an AdS'-Schwarzschild black hole in six-dimensions [48] . 
This is expected since, as we have shown, this specific Einstein-dilaton system in 5-dimensions 
was simply a dimensional reduction of a 6-dimensional theory containing only gravity and a 
cosmological constant. In fact, the dilaton solution presented may be explicitly uplifted to 
an AdS'-Schwarzschild black hole in six-dimensions [66]. While we have focused on a circle 
reduction from six-dimensions, similar statements can be made for generic n-torus reductions to 
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